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ABSTRACT 

As part of our survey of galactic stellar hales, we investigate the structure and stellar populations of 
the northern outer part of the stellar halo in NGC 55, a member galaxy of the Sculptor Group, using 
deep and wide-field V- and J-band images taken with Subaru/Suprime-Cam. Based on the analysis 
of the color-magnitude diagrams (CMDs) for red-giant-branch (RGB) stars, we derive a tip of RGB 
(TRGB)-based distance modulus to the galaxy of (m-M)o = 26.58±0.11(d = 2.1±0.1Mpc). From the 
stellar density maps, we detect the asymmetrically disturbed, thick disk structure and two metal-poor 
overdense substructures in the north region of NGC 55, which may correspond to merger remnants 
associated with hierarchical formation of NGC 55's halo. In addition, we identify a diffuse metal-poor 
halo extended out to at least z ~ 16 kpc from the galactic plane. The surface-brightness profiles 
toward the 2-direction perpendicular to the galactic plane suggest that the stellar density distribution 
in the northern outer part of NGC 55 is described by a locally isothermal disk at z < 6 kpc and a 
likely diffuse metal-poor halo with F-band surface brightness of fiy > 32 mag arcsec"^, where old 
RGB stars dominate. We derive the metallicity distributions (MDs) of these structures on the basis 
of the photometric comparison of RGB stars with the theoretical stellar evolutionary models. The 
MDs of the thick disk structures show the peak and mean metallicity of [Fe/H]poak ~ —1.4 and 
[Fe/H]niean '^ ~l-7, respectively, while the outer substructures show more metal-poor features than 
the thick disk structure. Combined with the current results with our previous study for M31's halo, 
we discuss the possible difference in the formation process of stellar halos among different Hubble 
types. 
Subject headings: galaxies: individual (NGC 55) — galaxies: halos — galaxies: structure 
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1. INTRODUCTION 

Our understanding of how a stellar halo surrounding 
a disk galaxy like the Milky Way has formed is still 
enigmatic, because the observational information of this 
faint galactic component is yet limited, except for Lo- 
cal Group galaxies. Halo stars in the Milky Way are 
characterized by low metal abundance and high veloc- 
ity dispersion and this extreme nature reflects the early 
chemo-dynamical evolution of the Milky Way. Exten- 
sive analyses of these stars have revealed, e.g., the dual 
nature of the halo structure (inner flattened halo in pro- 
grade rotation and outer sp herica l halo in retrograde ro- 
tation) fe.g.. ICaroUo et aLll2OO70 . Also, recent growing 
observational evidence suggests that the Milky Way halo 
has formed from an assembly process of many subsys- 
tems, as deduced from stream-like halo substructures 
(e.g., Yannv ct al. 2000). Similar substructures have 
also been identi fied in MSl ' s hal o, including the Giant 
Stellar Stream (jibata et al.l 120071 ). Indeed, MSl's halo 
provides a clear external view of ancient stellar popu- 
lations: the halo is found to be e xtended more than 
150 kpc (jGuhathakurta et al.l 120051 ) and has character- 
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istic spatial s tructure and metallicit y distributions (e.g ., 
iKalirai et aO2006l: llbata et al.ll2007tlTanaka et al.l201CII) . 
Another Local Group galaxy, M33, seems to have a 
rather smooth halo, although its disk orientation pre- 
vents the more refined picture of the halo. Thus, de- 
tailed studies of old stellar halos in various-type galaxies 
provide us with important clues to the understanding of 
galaxy formation. 

Exploring stellar halos even beyond the Local Group 
is of great importance to accomplish our understanding 
of their generic nature and past formation history in dif- 
ferent disk galaxie s at different envir onments. Model in- 
vestigations (e.g.. lKauffmannnil996[ ) predict that each 
disk galaxy has been developed through a different for- 
mation and evolutionary path: the collapse epoch, star 
formation history, and assembly rate of subsystems are 
not common and each stellar halo is expected to hold 
a different morphology, as other components (bulge and 
disk) differ along the Hubble sequence. 

In this regard, an edge-on spiral SB(s)m galaxy, 
NGC 55, which is a member of the nearby Scu lptor group 
that consists of approx imately 30 galaxies ()C6te et al.l 
Il997t iJerjen et al.ll2000( ) , is an excellent test-bed for this 
study: it provides an external perspective of a nearby 
late-type disk galaxy and yet is close enough to resolve 
individual stars. For example, NGC 55 shows asymmet- 
ric extra-pl a,nar morpholog:y as derived from emission- 
line images (iFerguson et al.iri996[ ) and its nuclear source 
shows, based on far-infrared measurements, energetic 
star formation rate, which is of comparable luminosity to 
the bright star fo rmation regions at the ce nter of M33 and 
the Milky Way (jEngelbracht et al.ll200l) . Based on the 
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obser vational s t udy o f asymptotic giant branch (AGB) 
stars, iDavidgd (|2005[ ) found that there are no evidence 
for a young or intermediate-age component in the extra- 
planar region such as 2 kpc off of the disk plane, sug- 
gesting that the po pulation has ag e s of about 10 Gyr. 
On the other hand, iTikhonov et alJ (|2005i) found an ex- 
ponential spatial distribution of RGB stars between 2 
and 7 kpc from th e center alo ng the minor axis, us- 
ing HST/WFPC2. ISeth et all (2005a,b), as part of a 
snapshot survey of 16 nearby edge-on late- type galax- 
ies, observed the thin disk of'NGC 55 using HST/ACS. 
They determined the fundamental properties of the re- 
solved stellar population such as the distance, MD and 
spatial distributions of main-sequence, AGB and RGB 
stars in the inner part at z < 1.5 kpc. However, since 
all the previous studies using HST are restricted to a 
small ficld-of-view, the global structure of an outer part 
of NGC 55 is yet unknown. Therefore, we conduct a high- 
resolution, deep wide-field survey of the halo of NGC 55 
using Subaru/Suprime-Cam, to elucidate its fundamen- 
tal properties for the first time. 

The layout of this paper is as follows. In Section [2l 
we present observations and summarize procedures for 
performing photometry for our Suprime-Cam data. Sec- 
tions [3l ID and [5] are devoted to our results for quantita- 
tive analysis of the CMDs of NGC 55's stellar popula- 
tions, including the distance estimate using the TRGB 
stars, the detection of substructures in NGC 55's halo 
from the stellar population maps, the radial profiles of 
the resolved stars, and their MDs using the theoretical 
isochrones. In Section |51 comparing with previous stud- 
ies on stellar halos, we discuss the implications of our 
results and conclusions are drawn. 

2. DATA 

In this study, we use the Suprime-Cam imager 
(jMivazaki et al.l [2002t ) on the 8.2-m Subaru Telescope 
on Mauna Kea in Hawaii. Suprime-Cam consists of ten 
2048 X 4096 CCDs with a scale of 0'.'202 per pixel and cov- 
ers a total ficld-of-view of 34' x 27'. We have observed the 
north part of a stellar halo in NGC 55 in 2009 December. 
Figure [T] shows the location of our Suprime-Cam pointing 
containing the edge of NGC 55's thick disk suggested by 
ITikhonov et al.l (|2005f ) and the outer halo out to around 
17 kpc perpendicular to the galactic plane. In this study, 
wc define that a;-axis corresponds to the galactic major 
axis, and z-axis which is converted based on the incli- 
nation (80°) in a direction perpendicular to the galactic 
mid-plane. The observations were made with Johnson V- 
band and Cousins /-band filters. Exposure times of our 
targeted fields are 960 sec and 1800 sec in V and /-band, 
respectively. The weather condition was photometric but 
slightly poor, with seeing of around I'.'O. To estimate the 
number of foreground stars of the Milky Way and un- 
resolved background galaxies in each NGC 55's field, we 
also obtained the imaging data for the control field, which 
is located at the same galactic latitude (b ^ —76°) as the 
object field, but is about 3° away from it. 

We have performed reduction, photometry and artifi- 
cia l star experim e nts fol lowing the same manner as given 
in I Tanaka et al.l (|2010( ). The raw data were reduced 
in the standard procedures, with the software package 
SDFRED, a useful pipelin e developed to optimally deal 
with Suprime-Cam images (jYagi et al.ll2002t lOuchi et al] 



120041) . We then conducted PSF- fitting photomet ry using 
thelR AF version of the DAOPHOT-K software (|Stetsonl 
|1987() . Comparing input artificial stars with output ones, 
we evaluated incompleteness due to low S/N ratio and se- 
lection criteria based on DAOPHOT parameters. Then, 
the 50% (80%) completeness of the object field for each 
fiher is reached at Vjl^css = 25.61 (Vff^css = 25.03) and 
/50,(^j.j. = 25.07 (/^°GC55 = 24.40) in the nearest field to 
the NGC 55 center, while the completeness of the control 
field reached at VS^,^^, = 25.80 (T/cZtroi = 25.23) and 
^Sontroi = 25.27 (/«o„,„, = 24.61) . The typical mean 
magnitude errors at the 50% (80%) completeness limits, 
based on the simulated numerous stars, are of the order 
of cry - 0.3 (0.2) and ct/ - 0.3 (0.2). 

Brief information about NGC 55 and the references 
are given in Table [1] Reddening correction is applied to 
this Suprime-Cam fi eld based on the extinction maps of 
ISchlegelet all(|1998l ). and the lDean etaP (|1978! ) redden- 
ing law /;(!/-/) = l.ME{B-V) and Ai = 1ME{V-I). 
The central coordinate of NGC 55 is derived from 2MASS 
observations, in agreement with the peak position of the 
blue light and also the location o f maximum symmetr y 
of Hi velocity field as described in lHummel et al.l ()1986fl . 
In addition, the galaxy's circular velocity, inclination and 
position angle are derived from their neutral hydrogen 
and radio continuum observations. 
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Fig. 1. — The locations and field of view of our Subaru/Suprime- 
Cam field, overlaid on a red Digitized Sky Survey image of NGC 55 
covering about 60' X 60'. A rectangle shows an area in which 
we analyzed our Suprime-Cam data, corresponding to 0.22 square 
degrees. 



3. DISTANCE TO NGC 55 

FigureHJshows the typical CMDs for stellar-like sources 
in the inner and outer region of NGC 55's halo (left and 
middle) and for the control field (right), after removing 
extended sources such as background galaxies and cos- 
mic rays based on DAOPHOT parameters. The three 
CMDs have the same sky coverage of 0.022 square de- 
grees. The solid lines in the CMD s show theoretical 
RGB tracks from iVandenBerg et all (|2006f) for an age 
of 12 Gyr, [a/Fe]= -1-0.3, and metallicities (left to right) 
of [Fe/H]= -2.31, -1.71, -1.31, -1.14, -0.71 and -0.30. 
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TABLE 1 
Properties of NGC 55 



Property 


Value 


Reference 


RA (J2000) 


00'^14™53!6 


NED 


Dec (J2000) 


-39°ll'48'/0 


NED 


Morphological type 


SB{s)m:sp 


NED 


A/ 


0.024 


Schlcgcl ct al. (1998) 


Vcirc (kms 1) 


110 


Hummel ct al. (1986) 


Inclination 


80° 


Hummel ct al. (1986) 


Position Angle 


109° 


Hummel ct al. (1986) 


(m - M)o 


26.58 (2.1 Mpc) 


this work 



The dashed lines denote about 50% and 80% complete- 
ness limits as determined by artificial star experiments. 

The clearest feature in the inner region of NGC 55's 
halo is a metal-poor RGB with -2.3 < [Fe/H] < -1.1. 
Furthermore, at / ^ 22.5, we also recognize a tip of RGB 
(TRGB). The stars distributed in the brighter part than 
the TRGB magnitude are probably metal-rich and/or 
young stellar population in the thermally pulsating AGB 
phase (TPAGB). In contrast, it seems that the GMD of 
the outer region is rather similar to that of the control 
field, suggesting that there might be only a little intrinsic 
halo population in this outer field beyond z ~ 12 kpc 
from the NGC 55 center. 

The TRGB is a useful in dicator to estimate the dis- 
tance to resolved galaxies (jSalaris fc Cassisil |2005() . If 
there are a sufficient number of old and metal-poor RGB 
stars in a targeted field, the TRGB is easily detected as a 
sharp cutoff of the luminosity function (LF) with the ap- 
plica tion of an edge-detection algorithm, the Sob el filter 
(e.g., ISalaris fc Cassisill2005l : iTanaka et aLll2010[ ). 

Detection of the TRGB is shown in Figure [3l The 
LF around the TRGB magnitude is derived from the 
selected RGB stars in a lower metallicity range than 
[Fe/II]= —0.7 to minimize the effect of empirical cal- 
ibration of the TRGB magnitude. The LF is gradu- 
ally standing up at Iq ~ 21.7 because of the increase of 
TPAGB stars and steeply rising up at /q = 22.48 ± 0.05. 
The TRGB magnitude suffers from statistical uncer- 
tainty of the order of about 0.05 mag. In estimating the 
extinction-corrected apparent magnitude of TRGB, sys- 
tematic errors also arise in association with zero point 
uncertainties, aperture corrections, photometric errors, 
smoothing of the LFs and extinction law; the resultant 
total error is evaluated by an rms of these errors. We fur- 
ther conduct 1,000 Monte Carlo realizations of the data, 
such that the brightness of each object is re-distributed 
in the form of a Gaussian distribution with standard de- 
viation equal to the photometric error. It is found that 
the resultant deviation of the TRGB magnitude is only 
as small as 0.03 mag and thus reasonably small. 

On the assumption that the absolute /o-band mag- 
nitude of the TRGB is M/^^ib ^ _4j ^ q ^ f^^. 

metal-poor TRGB sta rs, the distance modulus (e.g., 
ISalaris fc Cassisil 120051 ) is 26.58 ±0.11 corresponding to 
2.1 ± 0.1 Mpc. This value is consistent with previous 
studies based on the same TRGB meth od as this study 
(|Seth et al.ll2005al: [fikhbnov et al.ll2005l) . 

4. STELLAR POPULATION MAPS 



In this section, we investigate how the stellar popula- 
tions of NGC 55's halo are distributed spatially out to 
about 17 kpc from NGC 55's center. To do so, we draw 
stellar density maps corresponding to each stellar popula- 
tion, based on th e matched filter met hod explained in our 
previous study (jTanaka et al.l[2010l ). First, we visually 
separate the two major structures (thick disk and diffuse 
halo) and some substructures from the following spatial 
distributions of stellar population of NGC 55. Then, we 
quantitatively discriminate the difference of their popu- 
lations based on their CMDs and MDs. 

Figure [4] shows density maps of RGB stars in the 
north halo region of NGC 55, divided [Fe/H] into two 
non-overlapping ranges (left panel: [Fe/H]< —2.01, right 
panel: —2.01 <[Fe/H], see also Figure [2]). Then, we im- 
pose several criteria to get rid of unpleasant contamina- 
tions unrelated to the RGB stars from our photometric 
catalog: for example, appropriate color, magnitudes be- 
tween the TRGB and 80% completeness. In addition, to 
remove the remaining background and foreground con- 
taminations, we subtract the constant value estimated 
based on the data of the control field from each bin of 
the stellar maps. The background levels for RGB and 
AGB stars are /^yback "^ >^0.8 ± 0.1 mag arcsec"^ and 
Myback ^ ^'^•S i 0-1 mag arcsec"^, respectively. 

4.L Inner Structures 

Both panels of Figure |4] show that NGC 55 has an 
asymmetric stellar structure in eastern and western re- 
gions, regardless of metal abundances. There is a some- 
what extended structure towards the z-direction in the 
eastern part. We nominally refer the eastern and western 
structures to as "Thick Disk East (TDK)" and "Thick 
Disk West (TDW)", respectively hereafter. TDE seems 
to be related to the disturbed eastern disk which is in the 
shape of a tadpole tail, as seen in Figure [1] This might 
suggest that an accreted satellite galaxy gave a structural 
influence oii the eastern part of the disk struc ture (e.g., 
lAbadi et all 120031: lYoachim fc Dalcantonll2006[) . 

In contrast, assuming that TDW, which is a high den- 
sity region in the north-western part of NGC 55, is undis- 
turbed, we regard this region as the normal thick disk. 
The scale height of the plausible thick disk of NGC 55 
is estimated as zo,thick ~ 1.6 kpc (see also Figure JH), 
which is con sistent with that estim ated by the previous 
HST study (jTikhonov et al.ll2005D . In fact, considering 
Kirc = 110 km s~^, the estimated scale height is in good 
agreement with the relation between galactic rotational 
velocity and scale-height of thin/thick disk constructed 
by Yoachim & Dalcanton (2006) (see their Fig. 9). 

Next, we compare stellar populations between TDE 
and TDW, based on their CMDs shown in Figure [5j 
There are well-populated RGB stars in both CMDs, sug- 
gesting that a dominant population in both structures 
is old and metal-poor with [Fe/H]< —1. Furthermore, 
it is remarkable that AGB stars exist above the TRGB 
magnitude of /q = 22.48 as also shown in the previous 
section. The AGB population of TDE seems to be more 
sharply rising at {V — I) ~ 1.6 compared to that of 
TDW. If this feature is real, it indicates that the stel- 
lar population of TDE is different from that of TDW: 
the former is possibly several Gyr younger than the lat- 
ter. However, color distributions of both TDE and TDW 
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Fig. 2. — [(V — /)o, lo] color-magnitude diagrams for stellar-like sour ces in the inner and outer regions of NGC 55's halo (left and middle) 
and the control field (right). The solid lines are theoretical isochrones IIVandenBere et al.| [2(306^ of age 12 Gyr and [a/Fe]= -1-0.3 spanning 
the metaUicity range [Fo/H]= -2.31, -1.71, -1.31, -1.14, -0.71 and -0.20. The dashed lines denote the full ranges of the 50% and 80% 
completeness levels. 



at 21.00 < Iq < 22.48 have a similar shape and peak at 
{V — /)o ^ 1.7. Also, the ratio of the number of AGB 
stars at 21.0 < Iq < 22.48 and 0.6 < (V^ - /)o < 2.5 
to that of bright RGB stars at 22.48 < Iq < 23.0 and 
0.6 < {V - I)o < 2.0 is evaluated as about 30% for both 
TDE and TDW, where the effect of contaminations is 
limited within their Poisson noise. Thus, population dif- 
ference between TDE and TDW would not be significant 
from the current CMD analysis, so that that bulk popu- 
lations in TDE and TDW may originate from the same, 
pre-existing disk component. 

4.2. Outer Structure 

Despite metallicity cuts, the stellar density rapidly de- 
clines to a background level beyond z ~ 5 kpc, according 
to the stellar density maps of Figure SI However, metal- 
poor stars are diffusely distributed in the outer part of 
the halo. The left panel of Figure |6] shows a CMD in the 
outer part of NGC 55 which is enclosed by the largest 
rectangle in the metal-poor stellar density map of Fig- 



ure m whereas the right panel presents a CMD in the 
control field with the same area of a field-of-view as the 
left panel. Although both of the CMDs contain numer- 
ous background contaminations below /q ~ 23, the CMD 
of NGC 55's outer part seems to have a metal-poor pop- 
ulation with [Fe/H] < — 2 like the diffuse outer halo of 
M31 and the Milky Way. The existence of this metal- 
poor population in the outer halo can be verified by 
comparing both CMDs at bright and metal-poor RGB 
part of {V — I)o ~ 1.2 and Iq ^ 22.8 (corresponding to 
dashed boxes in the CMDs). In fact, the number of the 
stellar objects inside this box in the outer halo field of 
NGC 55 (nhaio = 69) is more than those in the control 
field (^control = 43), i.e., well beyond the Poisson noise. 
This fact that the diffuse halo extends to the outer part 
of NGC 55 is consisten t with the sugge stion of the pre- 
vious HST study (Tikh onov et al.ll2005[ ). However, since 
its signal-to- noise ratio is not high (of ^ 3), further ob- 
servations at other wavelengths are needed to reduce the 



A Subaru/ Suprime- Cam Survey of NGC55's Stellar Halo 




120 - 



-60 h 

J I I I I I I I I I 

22.5 

lo (mag) 

Fig. 3. — TRGB detection in NGC 55. The upper panel presents 
the smoothed LFs as a function of the incompleteness-corrected 
and background-subtracted /o-band magnitude zoomed in around 
the TRGB magnitude. The lower panel indicates the Sobel filter 
response to the LP. The vertical dashed line shows the derived 
TRGB magnitude. 



contamination of unresolved background galaxies and to 
confirm the presence of this halo component. 

In addition, two regions showing significant excess of 
stellar density beyond the level of the diffuse halo can be 
identified in the north-east field of NGC 55 (right panel 
of Figure |4]) , which we refer to as Substructure 1 and 2 
hereafter. These halo substructures are more clearly pre- 
sented in Figure[71 which shows a background-subtracted 
surface brightness profile in the same direction as the ma- 
jor axis at 5.25 < z < 7.25 kpc: the surface brightness 
is estimated using the resolved RGB stars selected as 
(Vb,/o) < (25.03,24.40), i.e., brighter magnitude than 
80% completeness limit (see also Section [6]). The error 
bars are estimated from the square root of the number 
counts including both stellar and background flux. The 
profile shows the presence of two overdense regions be- 
yond the Poisson noise at x '^ 2.5 kpc and a; ^ 8 kpc 
bounded by vertical gray dotted lines. The CMDs of 
both substructures (Figure |8]) show a similar distribu- 
tion of RGB stars including the presence of metal-poor 
stars with -2.3 < [Fe/H] < - 1.7, ahhough these sub- 
structures are spatially separated as shown in Figure [71 
The difference of stellar population between these sub- 
structures will be more quantitatively examined based 
on the analysis of their MDs in Section [S] 

4.3. Surface Brightness Profile 

To further investigate NGC 55's global structure, we 
investigate the distribution of surface brightness for the 
resolved RGB stars using those brighter than 80% com- 
pleteness limit. Figure [S] shows surface brightness profiles 
converting the summed- up flux counts of selected stars to 
the surface brightness in mag arcsec"^, and subtracting 
the surface brightness of the control field, for which re- 
maining foreground and background contaminati ons are 
removed based on the statistical method (Tanaka et al.l 
I2OIOI) . The error bars are estimated from the square 



root of the number counts including both stellar and 
background flux. Black and gray circles show the sur- 
face brightness profiles at < x < 4 kpc (where stars 
in TDW dominate) and — 5 < a; < — 1 kpc (where those 
in TDK dominate), respectively. It is noted that in the 
latter profile, there is a prominent overdense structure 
(corresponding to Sub 1) at z ^^ 6 kpc, as also mentioned 
in the previous subsection. 

It also follows from the figure that the thick disk 
component is extended at least up to about 5 kpc, as 
represented by isothermal disk models (dashed lines): 
the vertical profile of a dis k structure is defined as 
(jvan der Kruit fc Searlill981[ ): 



(1) 



I](z) oc sech — 



zq 



where S(z) is the surface brightness or density at a 
position z above the midplane and zq is the scale 
height. As discussed in Section 14.11 the scale height 
of TDW is 



^O.RGB 



— 1635 ± 30 pc which is a typical 

value for a late-type galaxy with Vdrc = HO km s~^ 
(jYoachim fc Dalcantonl [2C)0CT. However, on the assump- 
tion that NGC 55 has slower rotational velocity like 
60 < Vchc < 90 km s'^ (e.g., [Puchc ct all llQQlD . it has 
a somewhat thicker disk. Nonetheless, it is notable that 
the thick disk profile in the eastern part (TDK) has a 
significantly larger scale height, z^^q^ = 2203 ± 25 pc, 
than in the western part (TDW) , implying that the east- 
ern thick disk has been heated up by the interaction with 
an accreted dwarf galaxy in the Sculptor group. 

Regarding the surface brightness profile of the bright 
AGB stars (fitted with dotted lines in the figure), we 



obtain the scale heights of z. 



w 

O.AGB 



969 ± 62 pc and 



-^0 AGB ~ 1341 ± 60 pc, respectively, for the north- 
western and north-eastern regions of NGC 55. As 
reported by previous stud i es of NGC 55 and other 
galaxies (|Seth et al.l l2005bl : iTikhonov et al.l l2005| ). the 
intermediate-age AGB stars in spiral galaxies show 
steeper spatial gradients in their number density than 
the old RGB stars and thus are almost absent in the 
outskirts of galaxies. This result strongly suggests the 
presence of an older component with a larger s cale height 
in the o uter part of the disk as investigated in lSeth et al.l 
()2005bt ). Furthermore, the old RGB populations have 
a scale height similar to typical thick disk components, 
whereas the intermediate- age stellar population has a 
scale height very similar to typical thin disk components; 
th is is consistent with the more systematic consideration 
of lYoachim fc DalcantonI (|2006f) . Therefore, assuming 
our AGB population more clearly traces the thin disk 
of NGC 55, the eastern part of the thin disk in NGC 55 
may have been dynamically heated by merging satellites, 
forming a thicker disk as seen in the eastern part. In ad- 
dition, the fact that the scale height we calculated based 
on AGB is somewhat large r than that of the m ore inner 
eastern disk estimated by ISeth et al] ()2005b[ ) supports 
that the outer intermediate-age disk is more significantly 
expanded. Therefore, NGC 55's disk may be extended 
toward a large height at around 4 kpc from the galactic 
plane. 

Finally, we investigate the remote out er region of the 
diffuse halo of NGC 55 at z > 8 kpc. ITikhonov et~all 
(|2005f ) reported that there is a transitional point from 
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Fig. 4. — Log-scaled matched filter maps to a limiting magnitude of /q = 24.4 and Vq = 25.0. The left panel shows the distribution 
of red RGB stars chosen with (V - I) < 3.0 and -2.01 <[Fe/H], while the right panel is for blue RGB stars with 0.0 < {V - I) and 
[Fe/H]< —2.01. The resolution of this map is 0.15° X 0.15° pixels, smoothed with a Gaussian kernel over 3 pixels. 





(V - I)q (mag) 

Fig. 5. — [{V — /)o, ^o] color-magnitude diagrams for stellar-like sources in TDE (left) and TDW (right) regions of NGC 55. The solid 
lines are theoretical isochrones (VandenBerg ct al. 2006) of age 12 Gyr and [a/Fe]= +0.3 spanning the metallicity range [Fe/H]= —2.31, 
— 1.71, —1.31, —1.14 and —0.71. The dotted lines denote the full ranges of the 50% and 80% completeness levels. 



the disk to the halo at 6.5 kpc based on the deviation 
of their estimated density distribution from the starhght 
(see their Figure 12). However, considering our surface 
brightness profile shown in Figure [7] and [HI it is likely that 
the deviation is attributable to the overdense substruc- 
ture of Sub 2, taking into account their HST/WFPC2 
pointings on the minor axis of NGC 55. In fact, as dis- 
cussed in Section 14.21 the outer halo of NGC 55 is sub- 
stantially diffuse, and seems to have a nearly flat profile 
below the background level, My back ~ ^^-^ ±0.1 mag 
arcsec"^. Taking into account this faintness and spa- 
tial profile of the NGC 55's halo, its total mass may be 
smaller than that of M31, like halos in other late- type 
galaxies such as M33 and LMC. It is noted that the pres- 
ence of more substructures at a more remote field from 
the galaxy center (z ^ 12 kpc) may be likely if the halo 



of NGC 55 also originates from accretion of many dwarf 
galaxies, as suggested for bright spirals like M31 and the 
Milky Way (jBullock fc Johnstonll2005l ). However, to as- 
sess this picture, further observations of the outer halo 
of NGC 55 are necessary. 

5. METALLICITY DISTRIBUTIONS 

In this section, we discuss the difference of stellar 
populations between the above-mentioned substructures, 
based on the comparison for the MDs of the RGB stars. 
To construct MDs, w e adopt the Victoria - Regin a theoret- 
ical isochrones from IVandenBerg et all ()20Q6f ) (see also 
I Tanaka et al.|[20Tol) . In accordanc e with the in t erpola - 
tion and extrapolation scheme of iKalirai et al.l ()2006D , 
we calculate the metallicity for each star in the same 
segment of the CMD, assuming stellar population with 
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Fig. 6. — [{V — I)o , Jq] CMDs for stellar-like sources in the NGC 55's halo (left) and the control field (right). The solid lines are 
theoretical isochrones IIVandc nBcrg ct al. 2006) of age 12 Gyr and [a/Fe]= +0.3 spanning the metallicity range [Fe/H]= —2.31, —1.71, 
— 1.31, —1.14 and —0.71. The dotted linos denote the full ranges of the 50% and 80% completeness levels. The existence of metal-poor 
population in the halo field can be identified by comparing the number of RGB stars within the dashed boxes between these panels (see 
text for more details). 
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Fig. 7. — The surface brightness profile in the same direction 
as the major axis at 5.25 < z < 7.25 kpc. Two overdense regions 
bounded by vertical gray dotted lines are regarded as substructures. 
Sub 1 and Sub 2, which are also denoted in the stellar density map 
(right panel in Figure |4]l. 



[a/Fe]= -1-0.3, age of 12 Gyr and the same distance mod- 
ulus of 26.58 estimated in § [31 Furthermore, for secure 
determinations of MDs, we select the targeted RGB stars 
having 22.48 < Iq < 23.48 and 1.05 < (^ - /)o < 2.2. 
These selection criteria allow us to remove a number 
of contaminations, such as AGB stars, young stars and 
background galaxies. After performing the interpolation 
procedure for our photometric data, we subtract the MD 
of the control field from that of each object field in order 
to remove the effects of the remaining contaminations. 



Since the MD of the control field shows a nearly flat dis- 
tribution within the Poisson errors in the effective metal- 
licity range of —2.8 < [Fe/H] < — 0.7 (gray histograms 
in Figure [TU)) , this subtraction procedure does not affect 
the determination for the intrinsic shape of NGC 55's 
MDs. In addition, the fraction of the contaminations is 
estimated as roughly 30% in Sub 2 held with the smallest 
statistic, whereas it is no more than 6% in TDE with the 
highest statistic. 

Figure [10] shows the resultant contaminations- 
subtracted MDs for the four substructures we have newly 
identifled in this study. The vertical error bars denote a 
nominal uncertainty in each metallicity bin as derived 
from the Poisson errors. For reference, we also plot the 
MD of the control field (consisting of 69 objects) by a 
gray histogram, which clearly shows a nearly flat distri- 
bution within the Poisson errors. The vertical dotted 
lines correspond to the values of mean (black) and me- 
dian (gray) metallicity in each substructure, which are 
also summarized in Table [2] It follows that the MDs of 
TDE and TDW with high stellar density have a more 
metal-rich peak ([Fe/H] peak ~ —1.4) than those of Sub 1 
and 2 with low stellar density. All the MDs show a some- 
what broad distribution ranging from [Fe/H]~ —3 to 
~ — 1. For comparison, we also plot the M D for the Giant 
South ern Stream observed in M31's halo (jTanaka et al.1 
l2010( l (gray dashed histogram) , which is clearly different 
from the MDs in NGC 55. It is also worth comparing 
with the MDs o f NGC 55's thin di sk ([Fc/H]Jh_in^ - -1.0) 
as obtained bv ISeth et al.l ()2005b) (see their Figure 11): 
although the MDs of TDW and TDE show a similar 
shape to the MDs of the thin disk, the former are system- 
atically more metal-poor than the latter, even taking into 
account the likely effect of adopting the different theoret- 
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Fig. 8. — [{V — I)o, lo] CMDs for stellar-like sources in the two su bstructure regions, Sub 1 (left) and Sub 2 (right), of NGC 55's 
halo. The solid lines are theoretical isochrones (VandcnB crg et al.ll200^ ) of age 12 Gyr and [a/Fe]= +0.3 spanning the metallicity range 
[Fe/H]= -2.31, -1.71, -1.31, -1.14 and -0.71. The dotted linse denote the full ranges of the 50% and 80% completeness levels. 
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Fig. 9. — The surface brightness profile in the perpendicular 
direction to the galactic plane at a restricted x range, for the north- 
western (black) and north-eastern (gray) region of NGC 55. The 
long-dashed lines denote the model distributions of old RGB stars 
in the form of a locally isothermal disk model with a scale height of 
^W^j^g = 1635 ± 30 pc (black) and Zq^gb = 2203 ± 25 pc (gray), 
respectively. The dotted lines denote those of AGE stars with a 
scale height of zY'^q-q = 969 ±62 pc (black) and zJ^j^qb = 1341 ± 
60 pc (gray), respectively. The background levels for RGB and 
AGB stars are IJ-y^^y. ~ 30.8 ± 0.1 mag arcsec"^ and /^^"^^(.jj ~ 
33.6 ±0.1 mag arcsec"^, respectively. 



ical model of stell ar evolution and so mewhat large pho- 
tometric errors in ISeth et al.l ()2005b[ ). This may suggest 
some population difference between our detected thick 
disk and the inner thin disk in NGC 55. 

We also note that the MD in the individual region 
of NGC 55 shows a slightly different profile from each 



other: the MD of Sub f seems to show a higher fraction 
of metal-poor stars than that of Sub 2 and the MD of 
TDE has slightly a smaller fraction of metal-poor stars 
with [Fe/H] < - 2 than that of TDW. These differences 
in MDs may reflect different stellar population in each 
substructure. To quantify the difference in these MDs, 
we employ a two-sided Kolmogorov-Smirnov (KS) test 
for the cumulative form of MDs. The calculated KS 
probability for the null hypothesis that all of the four 
populations originate from the same stellar population is 
only less than 1%: this is a maximum probability for all 
of the KS tests. Thus, these MDs are statistically dif- 
ferent, thereby implying that stellar population in each 
substructure is different. However, we note that we can- 
not get rid of the degeneracies between metallicity and 
age when using isochrones to infer metallicity. 

6. DISCUSSION AND CONCLUDING REMARKS 

Our Suprime-Cam observation of NGC 55 has revealed 
that this late-type galaxy holds extended galactic struc- 
tures well beyond its bright disk component, namely 
TDE, TDW, Sub 1, Sub 2, and a diffuse halo. We discuss 
here the possible origin of these extended components in 
NGC 55. 

As already shown in Section 5, the stellar population 
in each of these structures seems to be statistically dif- 
ferent from each other, as deduced from the analysis of 
these MDs. However, taking into account the spatial dis- 
tribution of Sub 1 and Sub 2 as well as the difference of 
stellar populations between TDE and TDW, there is a 
possibility that both Sub 1 and Sub 2 may be associated 
with some merging event which also gives rise to disk 
thickening in its east part. If this conjecture is the case, 
then the MD of TDE should be a combination of the 
MD of TDW and that of cither Sub 1 or Sub 2, which is 
however unlikely as it follows from Figure 10. Thus, the 
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Fig. 10. — Metallicity distributions for the thick disk (upper panels) and substructures (lower panels). The number of stars used to 
derive these are shown in the upper right corner of each panel. Gray histograms correspond to the MD of the co ntrol field consisting of 69 
objects and gray dashed histograms correspond to the MD of M31's Giant Southern Stream obtained bv lTanaka, et al.i (i2010), which is not 
sensitive to the [Fe/H]> —0.7 metallicities. The vertical dotted lines show the mean metallicity (black) and median metallicity (gray) of 
each field. 



TABLE 2 

Fundamental properties of the substructures 



Name 



(mag arcsec"^) 



[Fe/HUe 
(dex) 



Standard Div. 

(dex) 



[Fe/H]^ed 
(dex) 



Quartilo Div. Error 
(dex) (dex) 



Thick Disk East 
Thick Disk West 
Substructure 1 
Substructure 2 



29.0 ±0.03 

28.7 ±0.02 

30.8 ±0.11 
31.2 ±0.17 



-1.68 
-1.76 
-1.99 

-1.84 



0.54 
0.58 

0.59 

0.42 



-1.57 
-1.67 
-2.13 
-1.72 



0.37 
0.46 
0.43 
0.30 



0.14 
0.14 
0.13 
0.13 



difference in stellar populations between TDE and TDW 
may be attributed to another merging event, which is not 
related to the origin of Sub 1 and Sub 2. 

Next, we compare the stellar populations of NGC 55's 
substructures with those of M31's halo. The MDs of 



NGC 55's substructures show a very different shape from 
those of M31's Giant Southern Stream, thereby sug- 
gesting that stellar population is quite different between 
NGC 55 and M31. Furthermore, the MDs of other sub- 
structures in M31 such as Stream C and D with low stel- 
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lar density ([Tanaka et al.l l2010l ) are also different from 
those of NGC 55. These results may suggest the differ- 
ence in the formation processes of stellar halos among 
different Hubble types of galaxies. 

Previous studies of M31's stellar halo show a posi- 
tive correlation between metal abundance and surface 
brightness for substructures observed in M31's halo 
(|Gilbert et al.ll2009l: [Tanaka et al.ll2Q10D . i.e., tidal debris 
with higher surface brightness tend to be more metal- 
rich. This correlation implies that tidal debris with 
higher surface brightness may originate from more lumi- 
nous and thus metal-rich dwarf satellites, and/or more 
recent encounters. For NGC 55's substructures, a similar 
correlation can be seen as shown in Figure lll[ although 
the trend is systematically shifted toward a more metal- 
poor end than that for M31's halo. Thus, NGC 55's 
stellar halo may have formed through merging of rel- 
atively metal-poor dwarf galaxies compared to MSl's 
stellar halo. Provided that the mass-metallicity rela- 
tion of dwarf galaxies in the Sculptor Group is the sani e 
as that observed in the Local Group ICote et ahl ()2000t ). 
NGC 55's halo may have originated from less massive 
dwarf galaxies than those for the formation of M31 's halo, 
thereby implying that NGC 55 itself is less massive than 
M31 at the current epoch. This seems to be in agreement 
with the properties of M33 as well: this late-type galaxy 
has a less massive halo than M31 and has some faint sub- 
structur es with metal-poor and low surface brightness in 
its halo (jMcConnachie et al.l[2010l ). 

The current study thus suggests that galaxies with dif- 
ferent Hubble types have different properties of stellar 
halos as characterized by different metallicity and sur- 
face brightness in their substructures. However stellar 
halos in earlier Hubble types have yet been unidenti- 
fied, simply because there are only few such early-type 
galaxies in the local volume so that current telescopes 
are able to resolve stars. On the other hand, a variety 
of substructures in the form of tidal streams have been 
found around more distant galaxies, where it is note- 
worthy that morphologies of tidal streams be tween early 
and l a te type galaxies are quite differ ent (e.g.. lPeng et al.l 
I2002t iMartinez-Delgado et al.l |2010[ ). Resolving these 
substructures into individual stars by, e.g.. Thirty Me- 
ter Telescope and/or Extremely Large Telescope, would 
provide further insights into the formation and evolution 
of stellar halos as a function of the Hubble sequence. 

As part of our survey of stellar halos in external galax- 
ies, we have observed the north part of the stellar halo 
in NGC 55, using deep and wide-field V- and /-band 
images taken with Subaru/Suprime-Cam. On the basis 
of the analysis of the CMDs compared with theoretical 
isochrones, we have obtained the following major results: 

1. We have found that the stellar populations at 
z > 5 kpc above the disk plane are dominated by 
old RGB stars with lower metallicity than M31 
and the Milky Way. We derive a TRGB-based 
distance modulus to the galaxy of (m — M)o — 
26.58 ± 0.11(d = 2.1 ± O.lMpc). 

2. Based on the density map of NGC 55's stellar pop- 
ulations, we have found asymmetric thick disk fea- 
tures in the north region of NGC 55, which we re- 
fer to as TDE for a disturbed thick disk in the 
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Fig. 11. — The median metallicity against surface brightness for 
the substructures of NGC 55 (filled black circles) and M31 (filled 
gray circles). Black crosses denote the disturbed and undisturbed 
thick disk structures of NGC 55. Metallicity uncer tainties (±0.2 
dex) are derived from the error analysis given in I Tanaka et al.1 
II2010I ). It follows that tidal debris with higher surface brightness 
tend to be more metal-rich in both galaxies and that NGC 55's 
trend is systematically shifted towards a more metal-poor end than 
that for MSl's halo. 



east and TDW for an undisturbed thick disk in 
the west. We have also identified two overdense 
substructures at z ^ 6.5 kpc, which we refer to 
as Substructure 1 and 2 in this work. These sub- 
structures may correspond to remnants of merging 
events of small galaxies associated with the forma- 
tion of NGC 55's halo. 

3. Detailed comparisons between the photometric 
data for halo regions of NGC 55 and that of the 
control field suggest the presence of a diffuse metal- 
poor halo extended out to at least z ~ 16 kpc. 
However, it is yet unclear to what extent this faint 
component is actually distributed. 

4. The stellar density distribution perpendicular to 
the galactic plane of NGC 55 is described by a 
locally isothermal disk at z < 6 kpc and a dif- 
fuse metal-poor halo with /iy si 32 mag arcsec^'^ 
at higher z, where old RGB stars dominate. 

5. From the photometric comparison of RGB stars 
with the theoretical stellar evolutionary model, 
we have obtained the MDs in the extraplanar re- 
gions. Both TDE and TDW show the peak, aver- 
age and median metalhcity of [Fe/H]poak ~ —1-4, 
[Fe/H],„oan ~ -1-7 and [Fe/H]med ~ -1-6, respec- 
tively. In contrast. Substructure 1 and 2 show more 
metal-poor features than the thick disk structures. 
Furthermore, the low KS probabilities indicate that 
all the MDs are statistically different, suggesting 
that the stellar population in each substructure 
may have a different origin. 
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